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The behavior of a two level atom in a half-cavity, i.e. a cavity with one mirror, is studied within the 
framework of a one dimensional model with respect to spontaneous decay and resonance fluorescence. 
The system under consideration corresponds to the setup of a recently performed experiment [J. 
Eschner et. ai, Nature 413, 495 (2001)] where the influence of a mirror on a fluorescing single 
atom was revealed. In the present work special attention is paid to a regime of large atom-mirror 
distances where intrinsic memory effects cannot be neglected anymore. This is done with the help 
of delay differential equations which contain, for small atom-mirror distances, the Markovian limit 
with effective level shifts and decay rates leading to the phenomenon of enhancement or inhibition 
of spontaneous decay. Several features are recovered beyond an effective Markovian treatment, 
appearing in experimental accessible quantities like intensity or emission spectra of the scattered 
light. 

PACS numbers: 42.50.Ct, 42.50.Md, 32.80.-t, 32.70.Jz 
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I. INTRODUCTION 

The change in the behavior of atoms when the struc- 
ture of the "surrounding" field differs from that of free 
space is treated so far in innumerable works and is essen- 
tially the basic topic of cavity QED jl], ||, [| . Effects like 
modified decay rates of atoms in cavities were visible in 
various measurements P, 0, H 0, H B|. In this case, the 
atoms couple irreversibly to a large number of field modes 
and the problem can be treated within the framework 
of perturbation theory. This regime is therefore known 
as "low-Q" or perturbative cavity QED [|l0). Another 
area of high importance is its counterpart, the physics 
of "high-Q" cavities, where the atoms interact strongly 
only with one (or few) field mode(s). In this context, 
recent experiments include, e.g., the observation of atom 
trajectories in cavities storing merely one photon [l2|| . 
Furthermore, among other things, effects caused by the 
spatial structure of a field mode in a cavity were demon- 
strated High-Q cavities serve also as a testing 
ground for fundamental quantum mechanical effects like 
entanglement or decoherence jl5[ 

Besides some considerations on the spontaneous decay 
of an excited two-level atom we will mainly focus in this 
paper on the problem of resonance fluorescence in a half- 
cavity, i.e. a cavity with one mirror, where we pay special 
attention to the position dependence of the atomic dy- 
namics. To this end we will particularly consider a phys- 
ical system which essentially coincides with the setup of 
an recently performed experiment jlTj. Here, the radia- 
tion which is emitted by a laser cooled ion stored in Paul 
trap is partly collimated by a lens and reflected back by 
a mirror to the atom. The intensity of the scattered light 
was measured as a function of the mirror position lead- 
ing to an oscillatory behavior of the photon counting rate, 
proving the existence of inhibited and enhanced sponta- 
neous emission effects. In this case, where the atom is 
relatively close to the mirror, the observed effects can in 
principal explained by introducing some effective modi- 



fied (position dependent) spontaneous emission rates and 
level shifts. This could be done since the time the light 
needs to bounce back and forth between the atom and 
the mirror could be set essentially to zero (Markovian 
limit). The situation is more complicated when the dis- 
tance between the atom and the mirror is large. 

In this paper we will, among other things, particularly 
consider this case and it turns out that the dynamics 
of the atom can be described generally in terms of non- 
Markovian (delay-differential) equations. As we will see, 
the distance between the atom and the mirror influences 
the atomic behavior essentially on two scales: On the 
one hand, on a large scale, i.e. whether the atom is lo- 
cated far away from the mirror or very close to the mirror. 
This scale can be measured essentially by a dimensionless 
quantity Tqt where Tq is given by the width of the field 
spectrum (in case of vanishing laser intensity it is sim- 
ply the atomic spontaneous emission rate) and the time r 
the light needs for a round trip between atom and mirror. 
On the other hand, the atomic behavior varies also if the 
distance is changed on the scale of an optical wavelength 
given by u>lt. For example in case of a small atom- mirror 
distance (Tqt <C 1) the equations of motions become ap- 
proximately Markovian and the well-known phenomenon 
of enhanced or inhibited spontaneous emission (depend- 
ing on lolt), can be recovered. Thus, it is possible to 
describe the system by introducing effective spontaneous 
emission rates and level shifts. In general, however, the 
retardation of the time argument in the equations of mo- 
tion cannot be neglected. We will not consider in this 
paper effects arising in the case of extremely small dis- 
tances, i.e. of the order of wavelengths or smaller, be- 
tween the atom and the mirror j|, [n], [ll| [lj], 20 . 



In connection with cavity QED, in the broadest sense, 
the above mentioned delay-differential equations ap- 
peared already in some works. These include for example 
the analytical treatments of Milonni et al. 0] in which a 
single excited quantum system coupled to an infinite set 
of equally spaced discrete levels was considered, a system 
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which reminds of an atom in a cavity but without tak- 
ing into account some position depending effects. The 
latter problem was discussed in the framework of a one- 
dimensional model by Cook and Milonni 0] in case of an 
excited atom in front of a partially reflecting wall (mod- 
eled as a collection of two-level atoms) and in a Fabry- 
Perot resonator. This treatment is closely related to our 
discussion of pure spontaneous decay since we will recover 
the same equation of motion. Another treatment of this 
problem was given by Feng and Ujihara p3| by using an 
appropriate set of mode functions in order to account for 
non-perfect mirror reflectivities. Dung and Ujihara J2J ] 
finally examined an atom in a three-dimensional Fabry- 
Perot resonator. Although a delay differential equation 
was not explicitly formulated, retardation effects in the 
interaction of two atoms were also discussed in [^5[ |2(| . A 
partly numerical examination of an atom inside a spheri- 
cal cavity was given by Parker and Stroud |27| . Numeri- 
cal examinations include furthermore the work of GieBen 
et q/.p8[ and Buzek et al. J29| both treating an atom in 
a (one-dimensional) cavity whereas the latter also inves- 
tigates the presence of material media. The mentioned 
works have in common that recurrences of the atomic 
population take place for large dimensions of the cavity 
which is due to one-photon wave packets bouncing back 
and forth between the cavity walls. However, there is al- 
ways only one excitation contained in the system making 
the problem accessible for analytical considerations. 

In case of a (near-) resonantly driven atom the dynam- 
ics of the system is more complicated since the number 
of excitations increase continuously. The scattered radi- 
ation will be reflected by the mirror and re-interact with 
its own source, the atom. This situation reminds us of a 
feedback problem [3^, |l| where mostly the assumption 
of a negligible feedback time delay is made. However, 
the situation of large atom-mirror distances would corre- 
spond to a non-Markovian feedback |32|, . 

This article is structured as follows. In Sec. II we will 
reconsider the problem of pure spontaneous decay of an 
excited atom in the presence of a mirror while in Sec. Ill 
a continuous laser excitation of the atom is incorporated 
to our examination. We will discuss several limits includ- 
ing low and higher laser intensities and small distances 
between the atom and the mirror. Finally, a summary is 
given in Sec. IV and cumbersome formulas and calcula- 



tions are moved to appendices. 



II. SPONTANEOUS EMISSION 

In this section we will investigate the spontaneous 
emission of an atom at rest in the presence of a mir- 
ror initially prepared in the excited state. We will derive 
a non-Markovian equation of motion for this system. Al- 
though this derivation is related to the calculations in |2|] 
it will be discussed here, not only to introduce our nota- 
tion but also to present new, entirely analytical results, 
also with respect to spectral properties of the emitted 



light. Furthermore it will turn out that some calculation 
methods can be transfered and some features of this sys- 
tem are recovered when we include laser excitation in our 
considerations. 



A. The model 

As already mentioned we examine an initially excited 
two-level atom with transition frequency loq in the pres- 
ence of a finite size mirror where the light emitted in 
a certain solid angle fraction e is reflected back to the 
atom. This is achieved by a lens which collimates the 
radiation before it is reflected [T^| . The remaining emis- 
sion is not affected by the mirror. Thus, it is reasonable 
to consider the coupling of the atom to two reservoirs 
(or "channels" ) consisting of one-dimensional fields with 
standing wave field modes and running wave field modes, 
respectively (see Fig. |l|), i.e. the Hamiltonian in rotating 
wave approximation reads 

H = H - d{E\{L)a- + <r+Ei(L)) 

- d(E$(0)a-+a+E 2 (0)) (1) 



with 



dk hujka\.a,k 



J dfc hw k b\b k , 



Ei(z) — i j dfc oik sin(fcz)afc, k > 0, 

E 2 (x) = i J dk /3 k e ikx b k , fceR, (2) 

and cok = \k\c. In contrast to a cavity, here, the mode 
density of the mirror channel is continuous since only 
one boundary condition has to be fulfilled. The opera- 
tors cr+, cr_ are the usual raising and lowering operators 
of a two-level system with upper level |e) and ground 
state \g), a+ = \e)(g\, a- = \g)(e\, and a\, b\, a k , b k are 
creation and annihilation operators of a photon in the 
fcth mode of the different environments. The dipole ma- 
trix element d is assumed to be real and for the sake of 
simplicity we suppress the vectorial character of d and 
E. The exact form of the factors a k and fit is of no 
importance here, we merely assume that they are ap- 
proximately constant in a frequency range of relevance 
(usually they have a frequency dependence a k ~ \p^k 
and (3 k ~ t/vk)- In order to investigate the dynamics of 
the system we make the Wigner-Weisskopf type ansatz 

m))= &e(*)|e,{0}i,{0}a) 



dfc&L(*)|ff,{*h.{0}2 



d*62.*(*)lff, {0}i,{A:} 2 



(3) 



where |{0}) denotes the vacuum state of the radiation 
field and |{fc}) the state with exactly one photon in mode 
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FIG. 1: Sketch of the physical system under consideration. 
The radiation emitted by a two-level atom is partly reflected 
(via a lens) back to the atom which is modeled by an atom 
coupled to two one-dimensional channels with a different 
mode structure (running and standing wave modes, respec- 
tively) . 



k. We will consider here an initially excited atom in the 
absence of any photon, i.e. b e (0) = 1 and b g k (0) = 0. In 
contrast to the notation in Eq. (|J), in the following the 
amplitudes are always taken in a rotating frame, i.e. we 
make the substitutions b e (t) — > b e (t)e~ luJot and b g k (t) — > 

b J g k (t)e~ iUkt . With the help of the "essential states" 
contained in the above ansatz it is possible to write down 
a closed set of equations of motion for the amplitudes 
which take the form 



dkg k sm(kL)e- i( -^° ]t bl k (t) 
dkh k e-^-^ni k {t), 



bl k (t)=g kS m(kL)e i ^-"°H e (t), 
i>lk(t) = h k e^~^%{t) 



(4a) 

(4b) 
(4c) 



with g k = a k d/h and h k = (3 k d/h. 

By formally integrating the latter two equations and 
inserting them into the first one we get 

b e {t) = -J*dt'b e (t') j dkglsirPikLy^-" o)(t '-* ) 

r dt'b e (t') f dkhle^-^W-V 

-7 / dtb e {t')5(t' -tfe-^oW-Q 



+ej / dt'b e (t')[S(t' -t + r) 
+8(t' - t - t)]c 



-iuJo(t' — t) 



(5) 



Here, we introduced the free space spontaneous decay 
rate 7 which is split up into a part £7 = irg ka /c and (1 — 
e)7 = 4irh ko / c, corresponding to the coupling of the atom 



to the first and second channel. The quantity e is the 
solid angle fraction which is covered by the lens since it 
characterizes the fraction of radiation which is reflected. 
In the above equations we also introduced the time r = 
2L/c the light needs for the distance atom- mirror-atom. 
Furthermore, in the first step of a Wigner-Weisskopf 
type approximation was made based on the well-known 
fact that the relative variation of g k , h\{~ u>) is very slow 
in the domain where the double integration in the first 
line of (||) lead to appreciable values. Diverging terms 
connected with level shifts are omitted. It should be 
stressed here that in this paper only distances between 
the atom and the mirror are considered which are much 
larger than an optical wavelength, i.e. loqt 3> 1 which is 
in an optical frequency domain already the case for, say, 
a millimeter. 

Eq. (||) yields finally the delay differential equation 



b e (t) 



r b e (t-T)0(t-T), (6) 



where Q(t) is the Heaviside step function. The first term 
on the right hand side of this equation corresponds to 
the usual free space exponential decay while the second 
term represents the effect of the reflected radiation on 
the atom which was emitted at time r before it interacts 
again with the atom. Thus, the retarded argument of the 
excited state amplitude directly indicates the memory 
effects which are inherent in the system. Furthermore, 
the second term is weighted with the factor e revealing 
that only a fraction of the emitted light is reflected. 

Eq. (|6|) is a delay- di fferent ia l equatio n [04| a nd since we 
will encounter in Sec. Ill B 2| and Sec. IIIB3 ; 



some more 



complicated equations of this type we move a further 
discussion of some general properties of equations of this 
kind to these sections. 



B. Discussion 

Using Laplace transformation and geometric series ex- 
pansion Eq. can easily be solved and one obtains 



be(t) = 



(£7/2)' 



n=0 



■e iu ° nT e-i^- nT \t - nT) n G(t - tit). 



(7) 



It should be mentioned that this expression can also 
be obtained by a direct Laplace transformation of the 
Schrodinger equation dJ) p5j . 

The above solution reveals that the systems dynamics 
has a "step" character which can be seen most easily if 
one divides the time axis into intervals of length r. For 
t £ [0, r] the sum consists only of one term, exp(— jt/2), 
which coincides with the free space behavior of a decaying 
atom. The physical reason for that is that the atom 
requires at least the time the light needs to get from the 
atom to the mirror and back to the atom again in order 
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to "see" the mirror. For t £ [r, 2r] the amplitude consists 
of two terms, 



b e (t) =c-^ + 



2 



T) (^-r), 



(8) 



giving rise to an interference term in the probability of 
finding the atom in the excited state. The second term is 
due to the emitted radiation reflected back to the atom. 
The light the atom emits right now arrives the atom again 
at the beginning of the third time interval where the sum 
in ^ includes a further term and so on. 

The role of the interference terms in the excited state 
probabilities strongly depends on the distance between 
the atom and the mirror which can be measured by the 
quantity jt. If we consider again the second time interval 
it is easy to see that that the interference term is of order 
7T while a further term is of order (jt) 2 which can be 
neglected for jt -C 1 (small distance). Hence we get the 
expression 



-T t (l+£7(t-T)cos(w T)), te [t,2t], 



(9) 



where we guess the beginning of an exponential series. 
The examination of the dynamics in this limit for larger 
times based on Eq.(Q) is relatively complicated. It is more 
convenient to return to the delay differential equation (^) . 
Since we are working in a rotating frame the amplitude 
b e (t) varies slowly on a time scale given by I/7. Thus, 
in the limit 77 C 1 we can make the approximation 
t — ► +0 in the argument of b e in the second term on 
the right hand side of Eq. (||) and obtain the Markovian 
equation 



b e (t) 



t < T 



.2(i_ £ e^ T )& e (i) ;t>r. 



(10) 



This leads to the excited state probability 



-It 

-7T p -7(*-t) 



t < T 
t > T 



(ID 



with 7 = 7(1 — £cos(o;ot)). The upper state population 
based on the exact amplitude (Q) in this limit is shown in 
Fig. ||. The behavior of the curves coincide almost per- 
fectly with the predictions of Eq. ( |i~l| ) : After a period of 
length r there is an enhancement or inhibition of sponta- 
neous decay depending on the factor 1 — e cos(wot) which 
corresponds to the amplitude of a standing wave mode 
sin (koz) at the position of the atom: In a node of the 
standing wave, spontaneous decay is inhibited while in 
an antinode it is enhanced. 

To get a more physical insight in this behavior the 
intensity of the electric field E\ which is reflected by 
the mirror depending on space and time is shown in 
Fig. ||. The derivation of an analytical expression for 
(E\(z, t)Ei(z, t)) can be found in Appendix ^ In a re- 
alistic situation with regard to the setup considered here 




FIG. 2: Upper state population of an atom close to the 
mirror for different exact positions, i.e. u)qt = 2nn (node) 
and ujqt = (2n + 1)% (antinode). The remaining parameters 
are e = 0.4, yr = 0.4. Also indicated is the corresponding 
free space solution (e = 0). 
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FIG. 3: Intensity of the field at different points of time for 
the parameters of Fig. ^| The atom is located at a node (i.e. 
loqt = 2n7r) at z = L. 



the frequency of the oscillations which can be seen in 
Fig. H (and also in Fig. ||) would be significantly higher 
than indicated in these figures. However, for the sake of 
visibility, a rather small frequency is chosen here. Due to 
the small distance between the atom (located at z = L) 
and the mirror (located at z = 0) the reflected light has 
the possibility to interfere with the radiation which is still 
emitted by the atom leading to a standing wave pattern 
of the form ~ sm 2 (koz) that has, in case of the example 
shown in Fig. 0, a node at the position of the atom, i.e. 
a zero electric field. Due to this fact, further emission of 
radiation in this channel is suppressed. Another interest- 
ing feature with respect to Fig. |^ is that the amplitude 
of the standing wave decreases for e < 1 with increasing 
time whereas the energy escapes in the other channel. 
The situation is reminiscent of a cavity where the atom 
acts like a partially transmitting mirror. 

In the limit of large distances between atom and mirror 
(i.e. 7T ^> 1) the sum (0) is dominated by the term with 
the highest power of jt. Hence, we get in a time interval 
[hit, (to + 1)t] 



\ to! 



3 -r(t-mT)( t _ mr jS 



(12) 



We see that the atom is partially reexcited by the ra- 



FIG. 4: Upper state population of an atom far away from 
the mirror and Wot = 2nn. The remaining parameters are 
e = 0.4, 7T = 10. The inset is a vertical magnification. 
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FIG. 5: Intensity of the field at different points of time for 
the parameters of Fig. ^j. The atom is located again at z — L 
while the inset is a vertical magnification. 



diation which it has emitted before and that the exact 
position (node or antinode) is not significant. This is il- 
lustrated in Fig. H and Fig. || where we plot again the ex- 
act solution for the excited state amplitude and the field 
intensity. The atom is placed in a node of the standing 
wave and the inset in Fig. [5] is a magnification in vertical 
direction of that part. We see that interference between 
outgoing and incoming radiation is much weaker than in 
Fig. | 

In this context it is also interesting to take a look at 
the spectrum of the emitted light (here this means the 
probability of finding a photon of frequency u> in the long 
time limit). This can easily be calculated by integrating 
Eq. @ and Eq. © and' using Eq. (Q) which leads to 



VJw,t) = 



AM 



% + i(u> - w) 



n=0 



I) 



■e iuJ71T (t - nr) n 



xG„h(i + i(wo - w))(t - nT)}Q(t - nr) 

(13) 

with 

G n [a] = iFi[n,n+l;s]-e*, (14) 
where iFi[n, m;x] is the confluent hypergeometric func- 



~o.s 



2.S 



FIG. 6: Transient photon population in channel 2 depending 
on frequency and time for a large distance between the atom 
(placed at an antinode, i.e. ujot = (2n + 1)tt) and the mirror. 
Further parameters are e = 0.4 and 7T = 10. 
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FIG. 7: Frequency dependent steady state photon popula- 
tion in case of a large distance between the atom and the mir- 
ror. The exact position of the atom is a node (ujot = 2n7r) 
or an antinode (loqt — (2n + 1)tt), respectively (e = 0.4, 
7 t = 10). 



tion and 



Aj(u) = 




f sin^r/2), j 
J = 2. 



(15) 



The transient photon population of the second channel 
in case of a relatively large atom-mirror distance for an 
atom placed in an antinode of the resonant standing wave 
field mode is shown in Fig. ^. In the long time limit the 
spectra take the form 



a 2 M 



ecos 



(lot)) 2 + (e^ sin(wr) + u — 



(16) 



Fig. shows the steady state photon population of the 
channel parallel to the mirror for different atomic po- 
sitions. The frequencies of the local minima in Fig. ^ 
and Fig. [7] coincide approximately with the frequencies 
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of those standing wave modes which have an antinodc 
at the position of the atom. This means that the pho- 
ton distribution in the second channel, which has initially 
roughly the shape of a Lorentzian (at the end of the time 
interval [0, r], see Fig. is affected by the back reflected 
light in the first channel. Here, the radiation components 
with the mentioned frequencies have a higher probability 
to be reabsorbed and emitted again (perhaps in modes 
of other frequencies) . This leads to a lower population of 
these modes. 

If the atom is very close to the mirror we recognize that 
the differential equation Eq. (|Io| ) contains the complex 
phase e luJ ° T , where the imaginary part of this factor can 
be interpreted as a level shift. This has consequences for 
the spectrum which takes in this limit the form 



The behavior of the system reminds us of that of a cas- 

, B7L pq, P9||, a formalism 



caded quantum system 35 



which deals with systems driven by non-classical types of 
light and which was applied in the theory of Markovian 



feedback 31 



The non-Markovian feedback contained in the system 
discussed here makes it difficult to solve the problem in 
an exact analytical way since it is not possible to estab- 
lish a closed set of equations describing the dynamics of 
the atom as in the theory of Markovian resonance fluores- 
cence. Thus we are restricted to approximative methods 
in the following sections. 



A. Perturbation theory 



l& 2 »| 5 



7 2 /4+^-w ) 2 ' 



(17) 



(18) 



where 7 is defined in Eq. (mh and 



u) = uj - e-sin(cj r). 



This expression can be derived with the help of Eq. ( ]1C 
(where we neglect the small contribution arising from 
the first time interval [0, r]) or with Eq. ( [l6| ) using the 
fact that for jt <C 1 the trigonometric functions in the 
denominator of Eq. ( |l6| ) vary very slowly on a frequency 
scale 7. The form of this spectrum illustrates first of all 
again the behavior shown in Fig. ||, i.e. the width of the 
Lorentzian is larger or smaller depending whether the 
atom is placed in an antinode or a node of a standing 
wave sin(fcoz). On the other hand the maximum of the 
function is shifted according to the imaginary part of 
the mentioned phase. The physical interpretation of this 
shift is based on the fact that the atom interacts with its 
own radiation. It corresponds to the energy of the atomic 
dipolc in the reflected electric field 
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III. LASER EXCITATION 

The system discussed in the previous section was 
amenable to an exact analytical treatment since the equa- 
tions of motion decoupled by using a Wigner-Weisskopf 
(or Markov)-type approximation. This reduced the prob- 
lem essentially to the solution of one equation describing 
only the atomic dynamics. This was possible because the 
system contained merely one excitation. The dynamics 
of the atom-field system becomes more complicated when 
we include a continuous laser excitation of the atom. The 
number of photons scattered by the atom into the two 
channels will permanently increase and a part of them 
will be reflected back, interacting again with the atom in 
addition to the laser light. The atom starts now to emit 
a different kind of radiation which again returns to the 
atom after some time and so on. Thus we expect that 
an electric field is constituted with a complex structure. 



Using the results of Sec. II we will, as a starting point, 
examine the influence of the laser for low intensities 
within the framework of a time dependent perturbation 
theory. Here, as well as in the following sections, the ef- 
fect of the laser is included in our considerations with the 
help of the standard semi-classical model for atom-laser 
interaction in rotating wave approximation, i.e. the new 
Hamiltonian reads 



H + V(t), 



with 



n 



V(t) = -h~-^-(e iuJLt a- 



-VUJLt, 



(19) 



(20) 



and laser- and Rabi-frequency u>l and f^o, respectively. 
Taking the ground state of the atom- field system as initial 
state we get in first order perturbation theory (assuming 
a weak laser intensity) the excited state amplitude in a 
rotating frame 



b L e {t) 



1 

ift Jo 



dt' (E\e 



-iH(t-t')/h v( ./\jHt/h 



i^-J dt'e- iAt 'b e (t'), 



V(t')e mt ' n \G) 



(21) 



with \E) = |e,{0}i,{0} 2 ), \G) = |s,{0}i,{0} 2 ) and laser 
detuning A = ujq — wl- The above expression is essen- 
tially equivalent to the one-photon amplitude of Sec. II B 



if the laser frequency ui^ is replaced by ui. Thus, we 
immediately get 



7 + 2iA 



n=0 



iG n [-a + iA)(t - nr)]9(i - nr). (22) 



Examples for the excited state amplitude are shown 
in Fig. g for different positions of the atom whereas the 
overall distance of the atom and the mirror is chosen 
to be quite large. The form of the curves has a direct 
interpretation. Since for low laser intensities coherent 
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FIG. 8: Upper state population in case of an atom far away 
from the mirror (yr — 20, e — 0.4) and weak laser excitation 
(fio = O.O57) f° r an atom placed in an antinode (dashed line, 
ujlt = (2n + 1)tt, A — 0), in a node (solid line, lult = 2nn, 
A = 0) and at a "slope" (dashed-dotted line, lolt = (2n—-^)n, 
A = -O.27). 



light scattering dominates the reflected radiation leads 
to a lower or higher "driving force" depending on the 
position of the atom. The system has similarities to an 
atom which is driven by two lasers where the phase dif- 
ference is controlled by the distance between atom and 
mirror. Actually the superimposed intensity of the lasers 
changes after each round trip of the light which leads to 
a transient upper state population as shown in Fig. ||. 
For example, if the atom is placed in a node the laser in- 
terferes always constructively with the "reflected" laser 
beam giving rise to a higher population in a time interval 
[nr, (n + 1)t] compared to the precedi ng one. This point 
will be further developed in the Sec. ( 1IIB 2 ) where we 
reconsider the limit of low laser intensities. 

The steady state population obtained from Eq. (|2l] ) is 
given by 



t^oo 7 £ _|_ 4A 2 

with modified decay rate and detuning 

7i = 7(1 - ecos(w L r), 
_7 



A 



A 



e-sin(o;Lr). 



(23) 



(24a) 
(24b) 



The situation is similar to the Markovian limit of the pre- 
vious section, i.e. we have a pronounced dependence of 
the atomic dynamics on the exact position of the atom 
(e.g. node or antinode of a standing wave of the laser 
frequency sin(fciz)). The difference is that this fact still 
holds in case of large atom-mirror distances. In the sense 
of Figs. (||) and (||) this is due to the fact that the in- 
terference ability of outgoing and reflected light does not 
depend on the distance since the laser provides a contin- 
uous scattered light field. 



A further quantity which is of interest in this context 
is the the second order intensity correlation function, 

(EjWEjit + T^it + T)^)) = \ aj \*G?\t,t + T), 

(25) 

where the index j = 1,2 indicates which channel is con- 
sidered and o.j = (5j2 — e)yh/(2d). Expressions for the 
electric field operators are given in Appendix |a|. This 
correlation function corresponds to the probability of de- 
tecting a photon at time t + T on condition that at time 
t a first one was detected (see e.g. pO|). 

Assuming that the detectors are arranged as in Fig. |ll| 
we have in the channel parallel to the mirror 



G 



{ 2 2) (t,t + T) 



(a+(jt)a+(t + T)a-(t + T)a-(t)) 
\\<T_U{t + T,t)<j-U{t,0)\G)\\ 2 . (26) 



By calculating the time evolution U in first order pertur- 
bation theory we obtain 



G¥\t,t + T) = \b^t)\ 2 \b^T)\ 2 , 



or in the long time limit 



lim G^\t,t- 



T) = - 



O 2 



7| 



4A- 



(27) 



(28) 



i.e. for 7T > 1 we have a behavior as shown in Fig. [g. 
This can be interpreted as follows. After the detection of 
the first photon the atom is in its ground state and has 
to be re-excited again before it is able to emit a second 
photon (antibunching). The radiation which is emitted 
now (the second photon) is split up into a part emitted 
into channel 2 and a part which is emitted into chan- 
nel 1. If T < r the light which is (or will be) reflected in 
channel 1 is not able to reach the atom before the photon 
detection in channel 2. Thus, we encounter, except for a 
constant factor, the same behavior as in free space. How- 
ever, if there is enough time for the radiation in channel 
1 to make a complete round trip, it is able to re-interact 
with the atom (in addition to the laser), which leads to 

a higher or lower emission probability in channel 2. 

(2) 

An expression for G\ ' (t, t + T) is derived in Ap- 
pendix W where we assumed again a detector arrange- 



ment as in Fig. 11, i.e. the atom is located in between 
the mirror and the detector. It is written as the norm of 
a sum of four states. These four contributions (in a non- 
rotating frame) are shown in Figs. ||a)-^|d) where each of 
these terms is connected to a different path leading to a 
coincidence detection at time t and t+T. The back action 
of the light on the atom is included in the dynamics of 
the time dependent operators. In principle the light has 
two possibilities to get to the detector: Either it takes the 
direct way or the indirect way via the mirror which leads 
to four possibilities for a two-photon detection amplitude 
indicated by the space-time diagrams in Fig. ^. For the 
sake of clarity the arbitrary distance between detector 
and atom, do, is not set to zero which is also indicated 
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«) 

a_{t+T-^)a.{t-^)\G) 



b) c) 



time 




time 




time 




d) 

-T^a.{t+T-^-r)a.{t-^)\G) 
time * 1 1 




FIG. 9: The four contributions to the second order intensity correlation function and the corresponding space-time diagrams. 
The correlation function is the squared norm of the sum of these vectors. The atom is located at L, the detector at zo and the 
distance between them is denoted as do- It is assumed that the detection of the first and second atom takes place at a time t 
and t + T, respectively. The »-symbol indicates the time of emission in the past relative to the detection times. The ordering 
of the operators in d) changes if T becomes smaller than r (dotted line) . 



in the time arguments of the operators. For calculations 
however, we will always set do = 0. The time arguments 
of the operators coincide with the emission time in the 
past relative to t and t + T, respectively. Note, that 
a possibility which includes a reflection gives rise to a 
negative sign and that the ordering of the operators in 
Fig. ^d) depends on the length of the delay interval T. 
The four contributions will interfere since they remain 
indistinguishable when a coincidence signal occurs. 

For weak laser intensities these quantities can be cal- 
culated using first order perturbation theory in a similar 
way as it was done in the derivation of Eq. ( p7| ) which 
leads to an expression of fourth order in Qq, 



Gf\t,t + T) 



\b L e (T)b L e (t)+e 2 ^b L e {T)b L e {t 
-e^b L e {T + r)b L e {t-r) 
-e^b L e {\T-r\)b^{t + S )\\ 



r) 



(29) 



where s = if T > r and s = T - r if T < r. We 
omitted step functions in this expression. In case of neg- 
ative arguments the corresponding quantities have to be 
set to zero. Recall that these amplitudes are written in 
a rotating frame. In the long time limit we have 



lim G^(t,t + T) = — 



O 2 
s-o 



4A 2 



x \2b L e {T) cos^r) - b L e (T + r) - b L e {\T - r|)| 2 . 

(30) 

This function is shown in Fig. [l(] for a relatively large 
atom-mirror distance. Several features are visible: First 

(2) 

of all G\ is not zero for T = 0. Indeed, diagram a) 
and diagram b) do not contribute to the detection prob- 
ability in this case which reflects the fact that after an 
emission process the atom is in its ground state and the 



+ 




FIG. 10: Second order intensity correlation function for a 
weakly driven atom (Oo = 0.057, e = 0.4) located at an 
antinode (lolt — (2n + l)7r), a "slope" (cult = (2n — 
and a node (lult = 2nn) of the standing wave sin(fcL2). The 
overall distance between the atom and the mirror is assumed 
to be quite large (7T = 20) while the laser is tuned to exact 
resonance (A = 0). 



probability amplitude that it immediately emits a second 
photon is zero. On the other hand, if the first photon is 
detected (and the atom is in the lower state) there is 
still the possibility that there is radiation around caused 

by a prior emission process which is represented by the 

(2) 

remaining diagrams. Due to this fact the value of G\ 
for T = differs from those for larger times where the 
"partial" antibunching effect of a) and b) decreases. For 
T = t we have a similar situation concerning diagram d) 
which does not contribute, i.e. a partial antibunching ef- 
fect which leads again to a different value of the detection 
probability compared to earlier or later times. 
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Within the framework of this perturbative treatment 
we can also calculate emission spectra without any effort 
which turn out to be monochromatic. However, since 
the results coincide with those of Sec. 
not quoted here. 



[II B 2 



they are 



B. Modified optical Bloch equations 

For further investigations it turns out that it is advan- 
tageous to work in the Heisenberg picture. 

As the amplitudes in the previous sections, in the fol- 
lowing the atomic operators and the mode operators 
are always represented in a rotating frame, i.e. <r_ 



-\u) L t 



er_ and a k , b k 



a k , e~ [uJkt b k . The Heisen- 



berg equations of motion for the operators a k (t), 

Ofc(t) = g k s\n{kL)<j_(t)e l ^"-"^\ (31) 
b k (t) = h k <r-(t)jl u *- u rt (32) 

yield after formally integrating and inserting into (g) and 
using a similar derivation as in (|5|) the electric field op- 
erator at the position of the atom, 

7 ih 



iuJLt (a-{t) - ee iuLT a-(t - r)6(t - t)) 
+ N 1 (t)+N 2 (t) (33) 



with E(t) = Ei(L,t) + (0, t) and noise operators 

' Pi r 

Ni(t) = jj &kg k sm{kL)a k {<S)e-™»\ 
N 2 (t) = ^ J dkh k b k {0)e-^ 1 . 
The second term on the right hand side of Eq. ( |33|) . 



-is^a-(t - r)e- lwi( *- T )6(t - r) 



(34) 



(35) 



can be identified as the reflected part of the electric 
(source-) field. Furthermore we will need in some calcu- 
lations the commutation relations of the noise operators 
and the "atomic" operators for t' < t, 



.7 



ih 



e-— e -^^-^S-i 
2d 11 



x [a^(t-T),a-(t')]e(t-T)G(t' -t + r). (36) 

Note, that for j = 1 this commutator is non- vanishing for 
f < t < t' + t (and t > t) in contrast to the Markovian 
case. With the help of expression (|3^), assuming again 
that all field modes are initially in the vacuum state and 
keeping normal ordering of the photon creation and an- 
nihilation operators, it is straightforward to derive a set 
of modified optical Bloch equations (OBEs), 



J 



At 
A 

d!< 

d 

dl ( 



<->=-(! + iA) (a-) i^(( C r +< 7_) - (a.a + )) - ((a + a.a.{t r)> - (a_a +( r_(i - r)))9(< - r), 



Ma 



i^((a + > - (<7_)) - 7 <<r+0 + el(e-^(a + (t r)a_> + e^(a + ^(t r)))@(t - r), 



dt 



dt {<7 - a+) = -dt {a+<7 - 



(37) 



r 



where we indicate for the sake of clarity only the retarded 
time arguments. As can be seen from these equations, 
the non-linear structure of the Heisenberg equations of 
motion leads to the appearance of correlation functions 
on the right hand side of Eq. (37). Hence, the modi- 
fied OBEs (|37]) can not be considered as a closed set of 
equations. However, it is convenient to take them as a 
starting point for approximative treatments. 



Small distance between atom and mirror (Markov limit) 



As in Sec. II B we will first consider the limit of small 
distances between the atom and the mirror, i.e. 7T <C 
1. Furthermore, we require now that the intensity and 
the detuning of the laser is not too high whi ch can be 
expressed by the condition fir <C 1 with il = ^/^q + A 2 . 



The latter defines a time scale on which the solution of 
the usual OBEsfe = 0) vary appreciably in the high 
intensity limit El), and we suppose that it will give us 
also an estimation of this scale in case of e 7^ 0. In 



Sec. Ill B 3 it will be shown that this condition has also 
some meaning in a frequency space. Thus, we can make 
again the approximation r — > +0 in the arguments of the 
operators in Eq. (5^) which leads to the equations 



df 
d 
dt 
d 
df 



= - If +iAW>-i^>, 



= - |l£-iA 



.n 0/ N 



(a z ) = iO ((<7+)-(^-))-7i(W + l), (38) 



where tl and A are defined by Eq. (p4|). For the sake 
of simplicity we set r = also in the arguments of the 



10 




| 




atom-mirror distance 



zo z 



FIG. 11: Sketch of the detector arrangement. A first de- 
tector opposite to the mirror (20 > L) measures the intensity 
7i in the first channel while a second detector measures the 
intensity I2 of the light in the channel parallel to the mirror. 



step functions (in contrast to Eq. (plf)), i.e. the differ- 
ence in the dynamics in the time interval [0, r] and later 
is neglected. Anyhow, taking the difference into account 
would merely lead to a (slightly) different initial condi- 
tion for Eq. (|3|) which would not alter the steady state 
results to be discussed here at all. 

Thus, the equations have the form of the familiar OBEs 
with modified spontaneous emission rate and detuning. 
They are rewritten as an inhomogeneous system of three 
differential equations where a z = <j + <j- — (T-<t+, since it 
is more convenient to perform steady state calculations 
in this representation. 

The steady state population of the upper state is easily 
obtained by inverting a 3 x 3-matrix which corresponds 
to Eq. ©, 



(0 + £7_ 



O 2 
"0 



7i 



4A 2 



Understood as a function of A this is essentially a 
Lorentzian (in the limit under consideration, treating the 
trigonometric functions as constants) with maximum at 
A max = sin(wLr) and width w — \/ 7 2 + 2 ffi . Apply- 



ing again the standing wave picture of Sec. fl h\ we see 
that the shift of the maximum vanishes if the atom is 
located at a node or an antinode of sin(cj^T/2). In con- 
trast to this, the width takes its minimum or maximum 
values at these points. Indeed we get a maximum shift 
if the atom is placed exactly in between a node and an 
antinode where the spontaneous emission rate is not al- 
tered at all. In the limit discussed here the steady state 
population is proportional to the measurable intensities 



ii ~ sin (wiT/2)(<T + ff- 

h ~ (o + <T_) ss 



(40a) 
(40b) 



of the light emitted in channel 1 or 2, respectively (see 
Fig. [ll] and Appendix |A|). Thus, a possible way to 




2a% 2(n+l)Jt 2(n+2)Jl 2(n+3)Jt 2(n+4)Jt 



FIG. 12: Excited state population depending on the distance 
between atom and mirror for A = (solid line) and A = 7/2 
(dashed line) and a larger solid angle (e = 0.4). The global 
distance is assumed to be small while the laser intensity is 
weak (fio = O.I7). 



demonstrate effects caused by the mirror is to measure 
the absorption spectrum of the atom, i.e. the intensity of 
the scattered light depending on the laser detuning. A 
further option would be the measurement of the inten- 
sity I2 for different positions of the mirror (keeping A 
constant), i.e. for different values of r, which would lead 
to a periodic variation in the measured intensity. This 
was done (for small global atom-mirror distances) in |l7j| 
whereas the measurement scheme slightly differed from 
that discussed here since the system under consideration 
was a three-level atom p7| . However, the basic principle 
is the same and for effects discussed in this paper it is 
sufficient to consider a two-level system. 



If we assume that e«l, Eq. 
lowest order in this parameter, 



I) can be expanded to 



(39) (<r+a-). 



fi 2 



l + 2e 



7 



/7 2 + 4A 2 



V 



cos(lult — (p) 



(41) 



with r = 7 2 + 2fl% + 4A 2 and tam(y) = 2 A/7. The 
presence of the relatively small level shift leads to a A- 
dependent phase shift (p with respect to the function 
cos(wlt) which corresponds to the phase of the standing 
wave sxsx(ujlt/2). The determination of this phase would, 
e.g., require the knowledge of the exact distance between 
atom and mirror. This difficulty could be avoided if one 
carries out a simultaneous measurement of I2 and I± since 
the phase of the latter is dominated for small e by the 
prefactor sin 2 (u> lt / 2) = (1 — cos(ujlt))/2. This means 
that, if there would be no level shift or A = 0, the two sig- 
nals were anticorrelated (i.e. a minimum of the Ji-signal 
coincides with a maximum of the /2-signal). The exis- 
tence of a level shift removes, in case of a finite detuning, 
this coincidence (see inset in Fig. [n]). For higher values 
of e we can get deviations from a pure sinusoidal behavior 
(see Fig. 
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2. Low laser intensity 



Provided with Eq. (|37]) we can now also examine the 
limit of small laser intensities in more detail. In this 
case (assuming that the atom is initially in the ground 
state) one expects that the atomic operators are approx- 
imately uncorrelated since coherent scattering processes 
dominate, i.e. we can make substitutions of the type 



(a q (t)a q ,(t'))n{a q (t)){a q ,(t')}. 
After some rearrangements one gets 



(42) 



dt 



|+iA) (cr_)-iin(i)<«r,), 



A ((T+(7 _) = i(n(t)<<7 + > -n*(t)<<7_» - 7 « - +ff _)) > 

(43) 

where we introduced the quantity 

H(t) = n - ie-fe^ LT (a-(t - r))0(i - r). (44) 

With the help of the decorrelation assumption (^2|) we 
eliminate the field degrees of freedom which are implic- 
itly still contained in Eq. ( |37j ) and get an equation for 
a reduced atomic system. This assumption is related to 
the fact that an atom initially in the ground state and 
weakly excited by a laser approximately behaves like a 
harmonic oscillator since a z = [(7+,er_] rs — 1. With re- 
gard to Eq. ( pt3| ) we have to replace (a z ) by —1 and it can 
be shown (see Appendix |b|) that assumption ([l||) holds 
in this case if the system is initially in the ground state. 

Anyhow, for the following discussion we keep for a 
short time the er 2 -term since the equations are more 
transparent in this form because the principal form of 
OBEs is conserved. Apparently the quantity H(t) can 
be interpreted as a modified Rabi-frequency in particu- 
lar if we recall the form of that part of the electric field 
operator which is due to the reflection of the light (see 
Eq. @), (E Tei (t)) « E (t)e~ iuJLt , with a slowly varying 
amplitude Eg(t). Thus, Eq. ( f44| ) can be written in the 
form n(i) = Qq + (2d/K)Eo(t), where the last term coin- 
cides with the definition of a Rabi-frequency. Let us as- 
sume now that jt 3> 1 and that the atom is in the ground 
state at t = 0. Then, the modified Rabi-frequency ( f44| ) 
has approximately the shape of a "stair function" (going 
up and down in general) with mostly decreasing distance 
between the single steps. This can be understood if we 
discuss the time evolution of the system in time intervals 
of length t: Between t = and t = r Eq. (^) are the or- 
dinary OBEs with Rabi-frequency Slo since the Heaviside 
function in (Q) vanishes. The solution of the ordinary 
OBEs yields for £1 < 7 



(M<)> 



i^o 



7 + 2iA ' 



te [0,r], 



(45) 



because the above expectation value is effectively con- 
stant after a few radiative lifetimes I/7 <C r. Accord- 



ing to this, in the next time interval [r, 2r] , the Rabi- 
frequency takes the form 



O = fio(l + fie [uJLT ) with /z; 



£7 



7 + 2iA 



(46) 



Now we have to solve again the ordinary OBEs which 
leads to 



(M*)> 



7 + 2iA ' 



te [r,2r], 



(47) 



giving rise to a new Rabi-frequency Slg in [2r, 3r] and so 
on. Thus, Eq. ( fi"3"l ) takes in every time interval [tit, (n + 
1)t], n E No the form of ordinary OBEs with different 

Rabi-frequencies fig™" 1 defined by 



T2 



(n-l) 



(48) 



while the expectation value of the dipole operator in an 
"intermediate" steady (cf. Fig. ||) state is given by 



n 



(n) 



7 + 2iA ' 



t G [tit, (n + 1)t]. 



(49) 



The value for n — * 00 of the recurrence relation (Q) is 
the limit of a geometric series or simply the fixed point 
of the map which is given by 



"0 — 



1 — fie 



1WLT ' 



(50) 



and leads to a steady state population (ct + (T_) ss w 

l^n° ' ) J 2 /(7 2 +4A 2 ) which coincides exactly with expres- 
sion ( 123 ) we got from the perturbation theory. 

Apart from this semi quantitative discussion it is 
worthwhile to study the equations of motion ( [43] ) in more 
detail. We explicitly make now the replacements cr_ — > c, 
0+ — » and therefore cr z — > —1, where c is a lowering 
operator of a harmonic oscillator. The equations of mo- 
tion then take the form 



| + .A)<e( t )>+^ 
+e l e i ^ T (c(t - r)>0(* - r), (51a) 



.n 



df 



<ct(t)c(i)> = -i-^« C (t))-( C t(t)))- 7 <ct(t) C (t)) 



(t-r)c(t)) 



C t(t)c(i-r)))e(t-r). 

(51b) 

In Appendix |b] it is shown that both (c^(t)c(t)) = 
(c 1 * (t)} (c(t)} and, as already mentioned, also the two-time 
correlation functions factorize if the atom is initially in 
the ground state. Furthermore, it is proved that these 
statements still hold in a steady state regime indepen- 
dently of the initial state. 
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Our task now is merely to solve Eq. (51a) assuming an 
initially unexcited atom. This equation is a linear delay 
differential equations, i.e. apart from the constant in- 
homogeneity, a type of equation like it already appeared 
in Sec. II. As a matter of fact, this equation reduces 
in every time interval [nr, (n + l)r] to a ordinary linear 
differential equation (ODE) with a time dependent in- 
homogeneity, so for any given initial state there exists 
a unique solution. In contrast to initial value problems 
concerned with ODE, delay differential equations need an 
initial function. In our case this initial function is defined 
due to the presence of the step function which yields in a 
first time interval [0, r] an ODE and is of of course given 
by its solution. This initial function is uniquely defined 
by the initial state and it replaces the quantity (c(t— r)) 
in the equation of motion in the next time interval [r, 2r] 
leading to an ODE with a time dependent inhomogene- 
ity As initial value we take of course the solution of the 
ODE in the first interval at t = r (which is justified since 
it can be easily shown that the solutions of the type of 
equations we consider here have to be continuous). The 
solution in [r, 2r] provides us again with the functions 
(c(i — t)) in [2r, 3r] and an initial value. Continuing this 
procedure, we see that we have to solve in every time 
interval [nr, (n + l)r] an initial value problem of ODE 
and we can apply all mathematical theorems which are 
concerned with such kind of equations. This "method 
of steps" J34| c an even yield analytical solutions as we 
will see in Sec. [II B J . Another method in order solve 



linear delay differential equations is by Laplace trans- 
formation like it was done in case of Eq. ( ]To| ) since in 
Laplace space the function with the retarded time argu- 
ment is simply replaced by the Laplace transformed of 
that function multiplied by an exponential function. 



Thus, the solution of Eq. (|43|) is unique (for a given ini- 
tial state) whereas the behavior of the derivatives is more 
complicated. With respect to this, it can be easily shown 
that (under certain conditions which are fulfilled in our 
case) the solution has at least nth continuous derivatives 
at t — nr and in general the (n + l)th derivative has a 
discontinuity. This feature can be identified in Fig 0, for 
instance, where we recognize a kink at t = r. 

In order to demonstrate the mentioned solution 
method we will now derive the transient solution of the 
delay differential equation ( |51a| ). The Laplace trans- 
formed of the expectation value (c) takes the form (as- 
suming (c(0)) = 0) 



£[(c(t))](iz + £) 



Q' 3 



1 



\z + £ iz + £ + at — a 2 e 



-(>*+Ot 



(52) 



£ G R+. We get 

/ m\ - ^ r a e(iz+i)t 1 

[C[t,) - 27rJ_ oo aZ iz + £ i_ a ? o-(*+«- iz + £ + Ql 



1 



Q'3 

2 



x / dz ■ 



^iz(t — nr) 



iz + £ (iz + £ + ai)" +1 ' 

(53) 

The Fourier transformation in the last line of this ex- 
pression (see e.g. E|) leads finally to the result fl22] ) 
already obtained from perturbation theory and a closer 
inspection of it confirms the result (|49|). 

The upper state population in a stationary regime co- 
incides therefore exactly with Eq. (^) which is the low 
intensity limit of Eq. (|39|). However, Eq. (j|^) is also valid 
for 77 > 1. The intensities in a stationary regime mea- 
sured in channel one and two, respectively (see Fig. |Tl| ) 
take again the form (^) which is due to the factorization 
property of the two-time correlation functions. Using this 
fact, we can furthermore easily calculate emission spectra 
of the light scattered in channel 1 or 2 which gives 

Si(w) ~ sin 2 (uj L T/2)\{c) ss \ 2 S(cj - cj l ), (54a) 
S 2 (iu) ~ \(c) ss \ 2 S(lu-lu l ). (54b) 

The fact that the spectra are monochromatic just ex- 
presses again that coherent, elastic scattering processes 
are involved in the limit of low laser intensities. 

It was already mentioned that the second order cor- 
relation functions factorize under certain circumstances 
but what about higher order correlation functions? A 
lack of the harmonic oscillator model is surely that in 
general the operator c(t) 2 is not equal to zero in con- 
trast to U-(t) 2 so we cannot necessarily expect that for 
example the quantity 

G ( 2 \t,t + T) = {c\t)c\t + T)c{t + T)c{t)) (55) 



with ot\ 



iA), a 2 



pla iuJ LT 
c 2 ' 



i^f-, and 



gives the correct result for T — > 0. In fact it can be 
seen from the results of Appendix [b| that G'' (t,t + T) 
is in general not equal to zero for T = 0, so as in the 
theory of ordinary resonance fluorescence (see, e.g., [[i"3"|) 
the correct result is obtained by perturbation theory. 

A further remarkable fact is that the harmonic oscilla- 
tor model reproduces the result of the Wigner-Weisskopf 
theory of Sec. II where pure spontaneous decay was con- 
sidered (see Appendix |b|) . At a first glance this seems to 
be surprising since the atom was initially in the excited 
state, i.e. (a z ) was far away from — 1. On the other hand 
we saw in Sec. II A that the state of the system is in this 
case always confined to the subspace spanned by the vec- 

tors {|e,{0}i,{0} a ),|ff,{fc}i,{0}a),|5,{0}i,{A;}a)} (and 
\g, {0}i, {0)2) if one wants to start in a state differ- 
ent from the excited state) which leads to the fact 
that the noise terms of the Heisenberg equations still 
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do not contribute to the modified OBEs. Further- 
more the two-time correlation function takes the form 
(a + (t')(T-(t)) = bt(t')b P .(t) and thus, the equation of mo- 
tion (from E g. (p7|) ) for the upper state probability is 



equal to Eq. (51b) for vanishing laser intensity. 



3. High laser intensity 

The examination of the systems dynamics for larger 
values of Qqt is more complicated since the incoherent 
nature of the scattered (and reflected) radiation becomes 
important. In order to investigate the dynamics in this 
parameter regime we will assume in the following that e 
is small so we can treat the "reflected" part of Eq. J37| ) 
as a perturbation. With the aim to obtain a closed set of 
equations which contain only terms of first order in e we 
can calculate the two-time correlation functions in 0-th 
order e depending on the initial state which is a single 
time expectation value and reinsert the result into (37). 



To this end we can multiply the Heisenberg equations of 
motion with cr+(t') from the left or U-(t') from the right 
where t' < t , make use of the commutation relations 
( |36| ) and calculate the expectation value. The equations 
we get in this way now contain third order correlation 
functions which are, however, of order e, and thus they 
are neglected. The solution for t' = t — r is given by 



C± (t, t - r) = U(t)C± (t-T,t-r) 



with 



d+(t,t') 



C-(t,t') 



(a + (t')a + (t)) 
(a_ («>_(*')) 

<<7 + (i)<7_(*')) 

V<<r_(*V+(*)M*')>> 



(56) 



(57) 



(58) 



The matrix elements Uij{r) of the evolution operator 
U(r) 



with 



A 4 = 





2 



iA 



ifio 



■Mo 
2 

Mo 
1 2 


h \ 

1 2 


-7 





7 


o J 



(59) 



are obtained by solving the corresponding differential 
equation. 

By inserting the 0-th order two-time correlation func- 
tions into Eq. (|37j ) we get finally an equation which is of 
first order in e, 

S(t) = A 4 S(t) + e K(r)S(t ~ t)9(< - r), (60) 

where we introduced the abbreviations 

S(t) ee «(,_(<)>, (a + (t)), (a + (t)a.(t)), (a_(t)<7 + (t))) T , 

(61) 



K(t) = 



with 





V i%/s(r) 







-i^/ 2 (r) 
if/*(r) 

7/4 (r) 
-7/*(r) 






0/ 



/iW = 
/ 2 (r) = 

Mr) = 

/ 4 (r) = 
Obviously, Eq. 



3 lW£T [C/34(r) 

.i^r 2i 7 



t/ 44 (r)], 



O 



^3l(T), 



e 1 — ^C/ 24 (r), 
"0 

^e" 1 — C/ n (r) 



(62) 

(63) 
(64) 

(65) 

^ t UI x {t)). (66) 



j) describes again a reduced atomic 

dyn amics b ut compared to the equation discussed in 
Sec. [Ill B 2 it is more complicated since we have now a 
coupled system of four delay differential equations. These 
equations are an extension of the ordinary OBEs (which 
are recovered in the e — > limit). We can apply the 
method of steps which yields the formal solution for times 
t G [mr, (m -I- l)r], 

S{t) = J7(t)5(0) + ^(e 7 ) m J7(t) / dti / dt 2 

n— 1 mr (m— l)r 

dt m B(t 1 )B(t 2 )...B{t m )S(0), 

(67) 

where £f(t) = U^ 1 (t)K(T)U(t - r). The above expres- 
sion has a form similar to that of the excited state am- 
plitude (fl0|). In fact, Eq. ( |67| ) yields in case of vanishing 
laser intensity 

(o- + (t)ff_(t)) = V^rfMtf-nr)" 

* * 71. ' 



71=0 



xe 



9(t - nr) (68) 

which is a acceptable approximation for e"/t -C 1. Fur- 
thermore, if we assume that jt -C 1, fir C 1 so that 
J7(r) w 1 and r — > +0 in the arguments of S, we recover 
Eq. (H) of Sec. ITIIBI . 



A numerically calculated example of the transient up- 
per state population is shown in Fig. [l^. 

The steady state solution can be found by calculating 
the eigenvector of the matrix A4 +eK(t) with eigenvalue 
0. From the form of the matrix K(t) we expect that 
the laser intensity (which is contained in the functions 
fj(r)) influences the decay rate(s) and driving force(s) in 
a steady state regime. In fact, we see in Fig. [l4] that the 
difference between the upper state po pulatio n obtained 
from Eq. (60) and the results of Sec. IIIB1 (indicated 
by the dashed lines) can be significant. Furthermore, for 
ilo 3> 7, small e and A = the upper state population 
takes approximately the form 



Q' 2 
"0 



1 + 2e^j- cos(w 7-)5( T ) 



(69) 
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FIG. 13: An example of the upper state population for fin = 
27, 7T = 5, e = 0.05 and A = 0. Until t — t the behavior 
equals to that of free space. In the magnified part (t > r) the 
functions are plotted for an atom in an antinode (lult = 2nir), 
a node (ujlt = (2n+l)7r) and for e = 0, i.e. no mirror (dashed 
line) . 
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FIG. 14: Excited state population deviation from the free 
space value (e = 0) in the long time limit depending on the 
distance between atom and mirror for A = 0, s = 0.1. The 
laser intensity is taken to be rather strong (fin = 2O7). The 
rapid oscillations have to be regarded in a rather symbolic 
way; in a realistic situation the frequency would be much 
larger. Also indicated are the maximum and minimum values 
(dashed lines) of the oscillations obtained by Eq. (fll|). 



with r = 7 2 + 20q. This expression equals Eq. (El]) 
obtained in the Markovian limit except for the function 
g(r) which is given by 



9(r) 



7 cos(f2 T) - — - sin(£l T 

4 Z7 



1 



We see that there is a modulation in the steady state pop- 
ulation defined by the Rabi frequency. This function has 
zero values in regimes fin'?" ~ Tin independently of uiqt. 
Thus, a strong laser can, in a way, inhibit the inhibition 
or enhancement of spontaneous decay. 

We will consider now the spectrum of the emitted light 
in the channel parallel to the mirror. For our purposes 
it turns out to be advantageous to define an emission 
spectrum in terms of the mean photon number increase 
N(t) of that channel in the long time limit, i.e. with the 



help of the differential equation (in a non-rotating frame) 

2, 



we obtain 



lim N(t) 

t — >oo 



with k,, 



hm I J dco (bt(t)bUt)) 
(l-e) 7 / dwSH, 



(71) 



where we defined the spectrum 



1 



S(u) = lim Re{(<7+ (*)&„(*)}}. 

nn tl , t— >oo 



(72) 



(73) 



The usual expression including the Fourier transformed 
of an atomic two-time correlation function is obtained 
(except for constant factors) by integrating Eq. (|7l| ) and 
inserting the result in Eq. (ff3|). Furthermore, corre- 
sponding to an operator O we define its fluctuating part 
60 = O — (O) with the help of which we can split the 
spectrum in a coherent and an incoherent component 



S(uj) = S coh (uj) + S inc (uj) 



with 



Scoh(w) 



1 



lim Re{(a + (t))(b u (t))}, (74a) 



lim Bje{(Sa + (t)5b u (t))}. (74b) 



It is easy to see that the coherent part of the spectrum 
takes the form 



Scon = (<7+)ss{V-)ssS(uJ - Wl), 



(75) 



where the stationary values are taken from Eq. (|60|). 

In order to calculate the incoherent component of the 
spectrum we can use a similar method as in the derivation 
of Eq. (|6^) . It is possible to derive a set of equations for 
the expectation values 

({Sa-(t)6b u (t)), (Sa+(t)6b u (t)), (6a z (t)8b„(t))) T = P(t) 

(76) 



(70) which takes in a rotating frame the form 



P{t) = [- i(w - u L )l + A 3 ]P(t) + K J (t) + ef(t,T). 

(77) 

Details of the calculation are given in Appendix [(]. The 
last term in the above equation includes two-time corre- 
lation functions which are again calculated in 0-th order 
e. This yields in the long time limit an expression of the 
form 



Ps. 



(78) 
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FIG. 15: Incoherent emission spectra for various distances 
between atom and mirror and weak laser intensity (Qo = O.27, 
A = 0, e — 0.15). The atom is always located in a node 
(ujlt = 2nn). 
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FIG. 16: Emission spectra for an atom for higher laser 
intensities (A = 0, e = 0.2). In Fig. a) the quantity 
Tor is small (fio = 37). The atom is placed in a node 
(solid line, ujlt — 2nir, 7T = 0.01) at a slope (dashed line, 
ujlt — (2n+i)7r, 7T = 0.005) and in an antinode (dotted line, 
ujlt — (2n + l)7r, 7r = 0.02). In Fig. b) the laser intensity 
is higher (£Iq = 57r7). Solid line: Node position, 7T = 0.1. 
Dashed line: Slope position, 7T = 0.1001. For visibility, this 
line is horizontally displaced by a small amount. Dotted line: 
Antinode position 7T = 0.1002. 



with 

M = -i(w - uj L )t + A 3 + ee- i( -"-"^ T K(T). (79) 

The atomic steady state expectation values which are 
contained in this expression are given by the steady state 
solution of the delay OBEs ( |30| ) . From this the spectrum 
Si nc can be calculated whereas for e = we get the usual 
Mollow-spectrum B ). 

Examples obtained from Eq. (fsh are shown in Fig. [l5| 
for weak laser intensity and an atom at a node. The 




FIG. 17: Illustration of the explanation for the different 
shapes of the emission spectra (see text) for increasing values 
of Tor. The curves in the figures correspond to the function 
sin 2 (u;T/2) for different values of r. For larger atom-mirror 
distances the relative variation of this function in a region To 
(the width of the fluorescence triplet) becomes significant and 
the Markov approximation is not valid anymore. The atom is 
placed at a node in this example. 



spectrum for 7T = in this figure is the Mollow result 
with a damping rate (1 — e)j. The structures arising 
at large distances resemble those of Fig. || and can be 
interpreted in a similar way. The situation changes in 
case of higher laser intensities. Examples are shown in 
Fig. 16 1) and |l6|b) for different values of r and different 
positions of the atom. We see that in general the width 
of the spectra varies and they are asymmetric depending 
on the position of the atom. 

This behavior can be understood at least on a qual- 
itative level if we take into account that a measure of 
the coupling strength of the atom to a field mode of fre- 
quency ui is given by sin 2 (wr/2). This function varies in 
frequency space on a scale 1/r. Defining the quantity 
To = 2f2 + 7 which approximately gives the overall width 
of the triplet we see that for T t <C 1 (and which we take 
now as the criterion for small atom-mirror distance) that 
the coupling is almost flat in the region where the spec- 
trum differs from zero (see Fig. p"7|a)). This situation cor- 
responds to the Markovian limit discussed in Sec. [II B 1 . 
Thus, we obtain in good approximation the usual Mol- 
low spectrum with a modified spontaneous emission rate 
7£. This is shown in Fig. |lqa) for various positions of 
the atom. The level shift ( |24bD , which acts here as a 
detuning in case of the dashed line, is so small that this 
curve cannot be distinguished from the Mollow spectrum 
with decay rate 7 on the scale of the figure. For larger 
values of Tor, but still 7T <C 1 , we have a situation like 
it is shown in Fig. |l7jb) where, as an example, an atom 
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located at a node of the standing wave sin 2 (wr/2) is cho- 
sen. For increasing laser intensity, the sidebands move 
towards regions of higher values of the coupling function 
leading to a higher damping of, say, the corresponding 
levels in a dressed state picture and thus to a broaden- 
ing of the sidebands (with increasing laser intensity until 
57 w 7t/t). For an atom placed at an antinode the behav- 
ior is simply the inverse. However, if the atom is placed at 
a "slope" , e.g. the one on the right hand side of the node 
which was considered in Fig. fl7|b), the spectrum becomes 
asymmetric since the transition responsible for the right 
sideband is stronger damped than the left one. Thus, the 
right sideband is broader than the left sideband which is 
in accordance with the dashed line in Fig. |l6|b) (For the 
sake of clarity, the dashed line in this figure is displaced 
horizontally by a small amount). The case Tqt 3> 1 is 
indicated in Fig. 17:) leading to structures as in Fig. [l5| 
or Fig. H 

So far we have discussed the case of exact resonance 
(A = 0) where the emission spectra are symmetric for an 
atom in a node or an antinode. This situation changes, 
in general, if we take a finite laser detuning. In case of 
Tqt <C 1 the spectra are approximately identical to the 
usual Mollow spectra with modified spontaneous emis- 
sion rate and detuning A, i.e. they are approximately 
symmetric independent of the exact atomic position. Ex- 
amples for this case are shown in Fig. |l8|a). Note that 
the sideband positions for an atom located at a slope are 
shifted towards the central peak which is due to the small 
frequency shift (the sideband positions are approximately 
given by ujl ± \/^o + (A — £7/2 sin(uj lt)) 2 ). This situ- 
ation differs from that when the distance between the 
atom and the mirror is increased. Here, the spectra be- 
come asymmetric even when the atom is located in a 
node or an antinode (see Fig. |T§|b)). 



IV. SUMMARY 

In this work we have discussed the behavior of an atom 
in the presence of a reflecting wall with regard to pure 
spontaneous emission, i.e. the decay of an initially ex- 
cited atom without any laser excitation, and with regard 
to an additional continuous driving laser field. In the 
first case, the one dimensional model applied here, can 
be solved exactly leading to a solution which directly 
reveals the retarded character of the system (photons 
bouncing back an forth between the atom and the mir- 
ror) visible in the state population, the field intensity 
and the (transient) photon spectrum. The limit of small 
distances yields the usual behavior of enhanced and in- 
hibited spontaneous emission which can be interpreted 
as an interference phenomenon of the the outgoing and 
reflected light pulse leading to a standing wave pattern in 
the field intensity: If the atom is placed in an antinode of 
this standing wave, spontaneous decay is enhanced while 
in a node it is suppressed. For large distances this inter- 
ference is not significant anymore and the node-antinode 
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FIG. 18: Incoherent emission spectra for non- vanishing laser 
detuning (A = IO7, e = 0.2, Slo = 37). Fig. a) shows an 
atom very close to the mirror: The solid line corresponds to 
node positions (cult — 2mr, 7T = 0.002), the dashed line to a 
slope position {lolt = (2n+i)7r, 7T = 0.0025) and the dotted 
line to an antinode position (ujlt = (2n + 1)71", 7T = 0.001). 
In Fig. b) Tor is larger. Solid line: ujlt = 2nn, 7T = 0.15, 
dashed line: ujlt = (2n + |)7T, 7T = 0.1505, dotted line: 
uj l t = (2n + l)n, "/t = 0.151. 



location of the atom becomes less important. The emit- 
ted photon wave packet is back reflected by the mirror 
leading to a partial re-excitation of the atom which starts 
now to emit again radiation and so on. 

In case of an additional driving laser the situation is 
more complex since the energy of the system increases 
continuously. Working in a Heisenberg picture we have 
derived a set of equations which serves as a starting point 
for several approximative treatments. In the limit of low 
laser intensities we saw with the help of perturbation 
theory and a harmonic oscillator model that the system 
behaves essentially like an atom driven by two monochro- 
matic lasers where the phase difference between the lasers 
is controlled by the atom-mirror distance. The intensity 
of the reflected light at the position of the atom depends 
on the intensity of the driving force on the atom at a 
preceding time which leads in general to a different state 
population in every time interval [mr, (m+l)r] (converg- 
ing to a steady state value). The dominance of coherent 
scattering was confirmed by the monochromatic emission 
spectrum of the system. In this limit we gave also a dis- 
cussion of the second order intensity correlation function 
which included in case of the field in channel 1 an inter- 
ference of different paths leading to a coincidence signal. 
This fact causes non-trivial structures in the correlation 
function. 

In case of a higher laser intensity incoherent scatter- 
ing becomes more significant. However, for small solid 
angles e it is possible to derive a closed set of linear de- 
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lay differential equations which represents an extension 
of the usual OBEs. It turned out that an intense laser 
field can significantly influence the system if we compare 
it with the Markovian limit where it is possible to de- 
scribe the system by OBEs with modified decay rate and 
transition frequency. With regard to the upper state pop- 
ulation, for instance, the laser can make the effect of the 
mirror to disappear regardless of the exact position of 
the atom (node or antinode). Furthermore, the influence 
of a strong laser was revealed by the emission spectra. 
Even if the width of the three peaks of the spectrum are 
each very small compared to the inverse delay time r, an 
intense laser field (or a high detuning) can "push" the 
sidebands of the triplet towards regions of a higher or 
lower coupling of the corresponding transitions to the ra- 
diation field leading to features like asymmetric spectra 
(see also @). 

A possible extension of the discussion presented in this 
article would be the inclusion of the motional degrees of 
freedom of the atom. Assuming an atom in a harmonic 
trap, as in the experimental realization fil} ], the reflected 
radiation will have an appreciable effect on the center 
of mass motion of the ion. This can serve as a further 
probe for effects discussed in this paper. Besides that, 
collective effects of two ions in the trap, like super- and 
sub-radiance, could be studied when the image of one ion 
is projected onto the other. The effect of one atom on an- 
other one, mediated by radiation over a large distance, is 
important for applications like quantum communication 
g§. 
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APPENDIX A: THE SCATTERED LIGHT FIELD 

Here we will sketch the derivation of an expression for 
the field intensity in the channel perpendicular to the 
mirror which was used for the generation of Fig. |j| and 
Fig. [|. Besides that, we will give formulae for the electric 
field operators in the Heisenberg picture which are used 
in the discussion of laser excitation and an expression for 
the second order correlation function which is needed to 
derive Eq. (|2^) . We use the coordinate system introduced 
in Fig. [I]. 

The intensity of the emitted light corresponding to 
Sec. I is defined by 

I(z,t) = (Ei(z,t)E(z,t)) = m)\EHz)E(z)m)) 



From Eq. (|) we see that 



A(z,t) = - e 



7 ih 1 
2~d 2vr 



dt'e-^ ot b e {t') 



x J dutpV-*) (e iw * -e- iuJ %) (c 

(A2) 

where the frequency integral gives rise to delta functions 
which yield non-vanishing terms only in certain regions 
of space and time, 

A(z, t) =siy^ 

x (e^i-iKit - § + §)0(t - § + f )0(f - §) 
+e -ia,o(f-§) be ^ + I _ i)e(i + | - §)©(§ - f ) 

_ e -o(t + J) 6e(< _ £ _ r )Q{t -i-^. (A3) 

The first and the second line of this expression repre- 
sent the outgoing light pulses to the right and the left 
side, respectively while the last line provides us with the 
reflected light pulse. 

In case of laser excitation the appropriate quantity in 
order to calculate the intensity is the electric field opera- 
tor in the Heisenberg picture. Starting with the Heisen- 
berg equations ( ^l| ) and the definitions (||) , we use a sim- 
ilar derivation as above, where A(z,t) has to be replaced 
by the operator Ei(z, t) and the amplitude b e by cr_ (and 
uo by u>l for a detuned laser). Apart from an additional 
noise term the result coincides with Eq. (A3). On con- 
dition that a photo detector is placed on the right hand 
side of the atom (cf. Fig. |ll|) at a position z > L the 
second line in Eq. ( |A3| ) does not contribute and one gets 



e 1ULT <7-{t- 
Nx(d ,t) 



r)0(* 



(A4) 



dkg k sm(k{d +L))a k (0)e- 1 ^ t , (A5) 



with 

Ni(d ,t) = - 
i 

where do = zq — L = zq — ct/2 is the distance between 
the detector and the atom. There are two different kinds 
of signals arriving at the detector, one which takes its 
way directly and one which takes the "loop way" over 
the mirror and therefore ne eds a lo nger time. 

If the conditions of Sec. Ill B 1 are fulfilled, we can 



approximately calculate the intensity in channel 1 by ne- 
glecting t in the arguments of the operators and the step 
function to obtain 



(E\{d ,t)E x {d Q ,t)) 



dfc a.k sin(fcz)6g k (t) 



= \A(z,t)\ 2 . (Al) 



(^) sinW2) 
x(a + (t-^)a-(t-^)) (A6) 
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which leads to expression (40a) if do is set to zero. 

For the sake of completeness we give here also the elec- 
tric field in channel 2 since it is used for various calcula- 
tions, 



E 2 (x,t) 



(l-e)7in 



-iui L t 



-(*--)©(*--)©(!) 



.(t+|)e(t+|)e(-£) )+N 2 (x,t) (A7) 



with 



N 2 (x,t) = -jdtW*(0)^ to - u,tl . (A8) 

With the help of Eq. ( |A4j ) one can easily find expres- 
sions for the intensity and the first order field correlation 
function in channel 1 (the functions connected with the 
channel 2 coincide with those of standard Markovian the- 
ory). 

Using Eq. ( A4) and the commutation relations (|36|) we 



get also an expression for the second order correlation 
function (p5|), 

G?\t,t + T)=\\ [<r_(* + ?>_(*) 

+o_{t + T -t)(t_(*-t) 

-e iWiT CT_(t + T)<T_(t-r) 



r) ff _(*)]|G) || 2 . 

(A9) 



We set the arbitrary distance do to zero and omitted the 
step functions in this expression which is valid if t > r (if 
not, components with negative arguments are simply set 
to zero) . The effect of the non vanishing commutator in 
Eq. ( 3q) is to conserve time ordering in the last term of 
Eq. (A9), i.e. the time argument of the operator on the 
left hand side is always greater than the right one. This 
is indicated by the symbol T_ . We see that for T < t 
and T > t the operators have to be exchanged. 



APPENDIX B: SOME FEATURES OF THE 
HARMONIC OSCILLATOR MODEL 

In order to derive an expression for the correlation 
functions in the harmonic oscillator model we start with 
the Heisenberg equation of motion for the operator c, 



dt 



c(t) = - a lC (t) + a 2 c(t - r)Q(t - r) + a 3 



id 



(Bl) 



with parameters on are defined in Sec. Ill B 2 and noise 



operators given by (p4|). Let us define a vector 

|*(t)> = c(tM0)>, i.e. |*(0)> = c(0)|V(0)), (B2) 



where IV'(O)) = \<Pi {0}i, {0)2) is the initial state of the 
system. The state 



\ ( p)=a\g) + b\e), |a| 2 + |6| 2 = l 



(B3) 



is an arbitrary state on the atomic space. Eq. (|Bl| ) pro- 
vides us with an equation of motion for this vector, 



-ai|¥(t)> + aa|¥(t-r)>e(t-T) 

+ a 3 |^(0)). (B4) 



Thus we have a linear inhomogeneous delay differential 
equation and its solution takes the form 



|*(i))=A(t)|^(Q))+B(t)|*(0)). 



(B5) 



This vector has the form |x(i))|{0}i, {0)2) where the non- 
constant part is an element of the atomic Hilbert space. 
The coefficients are given by 



A(t) =-V^(i- nr) n G n [- ai (t - nr)]e(t - r), 



B(t) 



n=0 



{t - nT) n e- ai{t - nT ^<d{t - nr). 



(B6) 
(B7) 



These expressions can be found by Laplace transforma- 
tion in a way it was demonstrated in Eq. (|5^) (the func- 
tion G n is defined in Eq. (|l4|)). 

Some expectation values of interest are 

(c(i)) = (V(0)|*(t)) = A(i) + B(t)<c(Q)>, (B8) 
( c t(i) C (t)> = |||*(t)>|| 2 

= (ct(t))(c(t)) + |B(i)| 2 (( c t(0) C (0)) - ( C t(0))(c(0)». 

(B9) 

From the above expressions it is immediately clear that 
{c\t)c{t)) = (J(t)){c(t)) Vt>0 (BIO) 

if (ct(O)c(O)) = (c t (0))(c(0)), which is the case iff the 
atom is initially in the ground state. In the long time 
limit this behavior is independent of the initial state, i.e. 



(c f c) 



(c) ss since limB(t) = 0. 

t — >oo 



If we start in the ground state the harmonic oscillator 
model reproduces result ( p2| ) gained from the perturba- 
tion theory and a remarkable fact is that Eq. (B9) also 
reproduces the solution (Q) obtained from the modified 
Wigner-Weisskopf theory if we set 03 = (no laser) so 
that A(f) = for all t and \<p) = |e). 

The two-time correlation functions take the form 

(ct(f)c(t)) = (*(f)l*(*)> 

= (ct(t'))(c(t)) +B*(OB(t)« C tc(0)) - (ct(0))(c(0))), 

(Bll) 
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where the time order is irrelevant. We see again that this 
quantity factorizes in special cases, 



J(t')c(t)) = (J(f))(c(t)), i£\ V ) = \g), 



lim (tf(t + T)c(t)) 

t — >oc 



= (ci) ss (c) SSl 



V|^). (B12) 



In order to derive fourth order correlation functions we 
proceed in a similar way. We define a vector 

|$(t0) = c{t')Mt)), (B13) 
|$(0)) = c(0)|*(t))=A(t)|*(0)). (B14) 



For this vector we get again an equation of the form (B4) 
where we replace |#(f)) -> |$(t')) an d |^(Q)) -> |*(t)> 
which finally gives the fourth order correlation function 

G { 2 2) (t,t + T) = {J(t)c\t + T)c(t + T)c(t)) 

= lll^ + T))|| 2 (B15) 

which leads to 

(t, t + T) = || (A(t)A(t + T) + A(t + T)B(t)c(0) 

+A(t)B(t + T)c(0))\^(0))\\ 2 , (B16) 

and thus, 

G^\t,t + T) = (c\t)c(t))(J(t + T)c(t + T)), \<p) = | ff ), 



lim Gf\t,t + T) 



(B17) 



From Eq. (A9) we also obtain 

.(2) 



\imGY>(t,t + T) = 16sin 4 (w L r/2)(c t c)i, V|y>). 

t—>oc 

(B18) 

This result is equal to the square of the intensity in the 
long time limit in this channel. 



APPENDIX C: CALCULATION OF THE 
SPECTRUM 

In this appendix a sketch of the derivation of the emis- 
sion spectrum in case of a higher laser intensity is out- 
lined. We will indicate in the following merely retarded 
time arguments. In order to get Eq. @ and Eq. @ 
we consider the operators 

8 A = 5<T-6b u = <T-b u + {cr-)(b u ) -bu(<r-) - a-{b u ), 
SB = 5a + Sb u = a + b 0J + {a + )(b 0J )-b u! (a + )-a + (b ul ), 
SC = 5a z Sb u! = a z b u + (a x ){b u ) -b u (a z ) - cr z {bj). 

(CI) 

After transforming in a rotating frame, 8A — > c~ 2lULt 8A : 
6B -> SB, SC -> e-' 1ULt 8C, ct_ -> e- iWit cr_, 6 W -> 
e~ lwt 6 w , the Heisenberg equations of motion for these op- 



erators yield, after taking the expectation value, Eq. (77), 



Pit) = [ - i(tu - Lu )t + A 3 ]P(t) + KuJ To(t) + ef(t,T), 

(C2) 








2 + iA 


V 


ifi 





where P(t) is defined in Eq. fl7q) and 



i . I 

-ifi 



I (t) = ( {a+<r-) - {a+)(<r-) 
-2<0+<7_><O 

The components of 7(t, r) are given by 



(C3) 
(C4) 



h(t,r) = - Ie u - ir ((5Ctr_(* - r)) + (tr,)ft_(t,r)) 



i"a(t,r)=-^e 



T ((a + (t-T)5C) + (a z )h + (t,r)) 

h{t,r) = 1 e- i ^ T ({a + {t-T)SA) + {a^)h + {t,T)) 
+ ^^((SBa^t - t)) + (a + )h-(t, r)) 



with 



h-(t,r) = (Sb u a-(t-T))^ UL - u ^ t , 
h+{t, t) = (a + (t - T)Sb u )^ UL -^*, 



(C5) 



(C6) 
(C7) 



where only retarded time arguments are indicated. The 
quantities h± can be calculated with the help of the 
Heisenberg equations of motion for b u . The result con- 
tains atomic two time correlation functions which have 
to be calculated again in 0-th order e depending on the 
initial state which is a single time expectation value. The 
calculation of the correlation functions contained in (C5) 
is analogous to the calculation done to derive the delay 
OBEs ( |60| ) : We multiplythe Heisenberg equations of mo- 
tion for the operators (CI) once from the left with cr + (t') 
and once from the right with <J-{t') (t' < t) and keep 
only terms of first order in e. The six equations we get in 
this way now contain again atomic two time correlation 
functions which have to be calculated as it was done for 
h±. Then we let t — ■> oo and get an expression of the 
form (Jt^) with 



K(t) = 



i/i(r) 







Vfr) 



-ifi / 3 (T) ifio/aW 



-i^/ 2 (r) 
7/4 (r) 



/iM = -e'^(r), 



Mr) = ^-e^g + (r), 



Mr) = -(e-^Ut^+e^U^r)), 



(C8) 

(C9) 
(CIO) 

(Cll) 

(C12) 
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and 



9 + {r) 



U(T) 





The inhomogeneity Ii(r) in Eq. 
we do not quote it here. 



|7q) is so lengthy that 



(C13) 
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